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1. It's all kinds of magic 

We offer a theory for counting magic squares and their innumerable relatives: semimagic and 
pandiagonal magic squares, magic cubes and hypercubes, magic graphs, and magical oddities like 
circles, spheres, and stars 

A magic square is an n x n array of distinct positive integers whose sum along any line (row, 
column, or main diagonal) is the same number, the magic sum. Magic squares date back to China 
in the first millenium B.C.E. [7j, came in the first millenium C.E. to the Islamic world and India 
[S], and passed to Europe in the later Middle Ages [S] and to sub-Saharan Africa not much after 
|17j . The contents of a magic square have varied with time and writer; usually they have been 
the first n 2 consecutive positive integers (standard squares), but often any arithmetic sequence 
and sometimes fairly arbitrary numbers. The fixed ideas are that they are integers, positive (or 
rarely, nonnegative) , and distinct. Even the mathematical treatises ^ |21 E] take positivity and 
distinctness so much for granted as never to mention them. 
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In the last century mathematicians took an interest in results about the number of squares 
with a fixed magic sum, but with simplifications. Thus diagonal sums were usually omitted and, 
most significantly, the fundamental requirement of distinctness was almost invariably neglected. 
Our work, however, follows tradition by adhering to the distinctness requirement. We shall call 
a square magic or strongly magic if its entries are distinct and weakly magic if they need not be 
distinct; and strongly or weakly semimagic if the diagonals are ignored. 

The literature on exact formulas examines the functions W(t) and Wo(t) that count positive and 
nonnegative weak semimagic or magic squares with magic sum t. These functions are amenable 
to analysis in terms of Ehrhart theory [5]. (See for example |1U1 113j for semimagic squares; for 
magic squares see 3 .) One treats a square with upper bound or magic sum t as an integer vector 
x £ t ■ [0, 1]" , confined to the subspace ts, where t = 1, 2, 3, . . . and 
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s := {x £ M. n : all line sums equal 1}, 

the magic subspace. (Exactly which subspace this is depends on whether we treat squares that 
are semimagic, magic, pandiagonal magic [where we include among the line sums the wrapped 
diagonals], or of another type.) Thus a square x is an integer lattice point in the t-dilate tP of the 
magic polytope P : = [0, 1]™ fls; moreover, x S tP°, the relative interior of P, if and only if the 
square is positive. 

Ehrhart's fundamental results on integer-point enumeration in polytopes Q] give much insight. 
Ehrhart theory implies that W and Wq are quasipolynomials in t. (A quasipolynomial is a function 
Q(t) = Y2o c rf l w ith coefficients q that are periodic functions of t, so that Q is a polynomial on 
each residue class modulo some integer, called the period; these polynomials are the constituents 
of Q.) The quasipolynomials are polynomials (that is, the periods are 1) in the semimagic case, 
because the matrix that defines s is totally unimodular so the vertices of P are all integral. In the 
magic case this is unfortunately not so and the period is not easy to calculate. 

Still there were no exact (theoretical) formulas for strong squares (not even in the comprehensive 
tome |12p. with the exception of Stanley's |15l Exercise 4.10]. With the theory of inside-out 
polytopes fl] we can attack this and related counting problems in a systematic way obtaining a 
general result about magic counting functions and an interpretation of reciprocity that leads to 
a new kind of question about permutations. In inside-out theory we supplement the polytope 
P = [0, l] n PI s with the pair-equality hyperplane arrangement 

•K := l K[K n2 ] s = {h i:j Ds:i<j< n 2 }, 

where hij is the hyperplane Xi = xj, "K[T] := {hij : ij € E} is the hyperplane arrangement of the 
graph r with edge set E, and denotes the complete graph on d nodes. The number of n x n 
squares corresponding to s with magic sum t is the number of integer points in t (P° \ |J !K). This is 
a quasipolynomial in t by the general theory of inside-out polytopes. Then inside-out reciprocity [I] 
gives the enumeration of weak nonnegative squares with multiplicity; this is reminiscent of Stanley's 
theorem on acylic orientations |14j . 

Another famous kind of square is latin squares and their relatives. Here each line has n different 
numbers. In a latin square these n numbers are the same in every line and are normally taken to 
be the first n positive integers. In any latin square in this broad meaning, every line has the same 
sum. Suppose we add this property to the definition of a latin square but we loosen the restriction 
on the entries, so that the square is filled with positive integers having equal row and column sums. 
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We call such squares magilatin. Then a magilatin square is a point in Z n ; the only difference 
between a semimagic and a magilatin strong square is that we assume fewer inequations between 
the entries; while in a semimagic square each entry must differ from every other, in a magilatin 
square it must differ only from those that are collinear with it, a line being a row or column. As 
with magic and semimagic squares, inside-out polytope theory yields theorems about the number 
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of magilatin squares as a function either of the magic sum or of the largest allowed value of an 
entry in the square. 

There is a parallel generalization of latin rectangles. A latin rectangle is an m x n rectangular 
array filled by n symbols, none repeated in a row or column. The asymptotic numbers of latin 
squares and rectangles of given dimensions have been the subject of many studies. Our geometric 
counting method leads in a different direction that, as far as we know, has not been studied. Define 
a magilatin rectangle to be a point in Z mn but not in |J J{[K m x K n ]. We discuss the number of 
magilatin squares or rectangles of fixed dimensions as we vary the maximum permitted value. 

The magic and latin properties generalize far beyond squares and rectangles. Semimagic and 
pandiagonal magic squares suggest a general picture: that of a covering clutter, consisting of a 
finite set X of points together with a family L of subsets, called lines for no particular reason, 
of which none contains any other and none is empty, and whose union is X. We want to assign 
positive integers to X so that all line sums are equal to a single number. Such a labelling is called 
a weakly or strongly magic or latin labelling of the covering clutter, depending on the particular 
requirements. There are many interesting examples that we cannot treat individually here. (X, £) 
may be a finite affine or projective geometry, the "lines" being the subspaces of any fixed dimension; 
more generally it may be a block design. It may be an n x n x • • • x n hypercubical array. It may 
be a k-net, where the lines fall into k parallel classes (with k > 2) so that each point belongs to 
a unique line in each parallel class. (A semimagic square is a 2-net and a pandiagonal square is a 
kind of 4- net.) All these examples have lines of equal size, a property that has advantages but is 
not necessary for the theory to apply. 

These ideas generalize still further. Take rational linear forms fx, fi, . . . , f m . A magic labelling 
of [d] := {1,2, ... ,d} with respect to fx, /2, • • • , fm is an integer point x £ M. d such that all the 
values fi(x) are equal to the same number. The analog here of a covering clutter in which all lines 
have the same size is a system of forms for which all values fi(l) (the weights; 1 is the vector of 
all ones) are equal; such systems have nice properties. 

We treat two distinct interesting approaches to enumeration. Traditionally, magic and semimagic 
squares have been counted with the magic sum as the parameter (due to its geometrical interpre- 
tation we call this affine counting); but another tack is to take as parameter the maximum allowed 
value of a label (which we call cubical counting) . These same two systems apply to latinity. In our 
treatment we develop both counting systems equally. 

The reader may wonder how practical our counting method is. We believe it is relatively feasible. 
In [Sj we apply it to solve in utter detail six problems of 3 x 3 squares: magic, semimagic, and 
magilatin (all strong), counted both cubically and affinely. 



2. Inside-out polytopes take the stage 

The theory of inside-out polytopes [I] is designed to count those points of the integral lattice 
Z d that are contained in a rational convex polytope P but not in a rational affine hyperplane 
arrangement "K, that is, where each hyperplane is spanned by the rational points it contains. We 
call {P, "K) a rational inside-out polytope, closed if P is closed. We shall always assume P is closed, 
except when we specifically state otherwise. The affine span aff P may be a proper subspace of R d . 

A region of IK is a connected component of M. d \ (J "K. A closed region is the closure of a region. 
A region of {P, "K) is the nonempty intersection of a region of IH with P. A vertex of {P, "K) is 
a vertex of any such region. Note that a closed region of {P, "K) is the closure of an open region 
of {P, "K) and therefore meets the relative interior P°. The denominator of {P,"K) is the smallest 
positive integer t for which t~ lr L d contains every vertex of (P, Jf). 
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The fundamental counting functions associated with (P, K) are the closed Ehrhart quasipolyno- 
mial, 



where P is closed and the multiplicity mp^c(x) of x E R rf with respect to IK is defined through 



mp^ix) :-- 



the number of closed regions of (P, "K) that contain x, if x € P, 
0, if x i P, 



and the open Ehrhart quasipolynomial, 

E°p« M (t) :=#(t- 1 z d n[P°\UK] 



We denote by volP the volume of P normalized with respect to Z d n aff P, that is, we take the 
volume of a fundamental domain of Z d D aff P to be 1. When P is full dimensional this is the 
ordinary volume. 

Theorem 2.1 (|IJ Theorem 4.1]). If (P,!K) is a closed, full- dimensional, rational inside-out poly- 
tope in M. d , then Ep^(t) and E° po ^Ai) are quasipolynomials in t that satisfy the reciprocity law 
Ep nr(t) = (—l) d Ep : ji(—t)., with period equal to a divisor of the denominator of (P, with 
leading term (vol P)t d , and with constant term Ep^(0) equal to the number of regions of (P, Ji). 

In particular, if (P, K) is integral then Ep^ and E po ^ are polynomials. 

The Mobius function of a finite partially ordered set S is the function [i : S x S —* Z defined 
recursively by 

if r ^ s, 

n(r, s) := < 1 if r = s, 

,-Er< M < s M r ^) ifr<s. 

Sources are, mfer alia, and ^3 Section 3.7]. 
The intersection semilattice of K is defined as 

L(K) : = {H§ : S C M and H§ ^ 0}, 

ordered by reverse inclusion (so the whole space M. d , the intersection of no hyperplanes, is the 
bottom element 0). Its members are the flats of "K. The intersection poset of (P°, K) is defined as 

£(P°, K) := {P° n OS : § C K} \ {0}, 

ordered by reverse inclusion. £(P° , IK) is a ranked poset and every interval is a geometric lattice. 
It equals £(IK) if f| 5£ meets P°. 

The arrangement induced by IK onsCl d is 

IK S := {ft n s : /i e M, ^s}. 

For the second theorem we need the notion of transversality. IK is transverse to P if every flat 
u € £(IK) that intersects P also intersects P°, and P does not lie in any of the hyperplanes of IK. 
Let 

E P (t) := # ftPnZ d ' 



the standard Ehrhart counting function (without any hyperplanes present). 
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Theorem 2.2 (|IJ Theorem 4.2]). If P and "K are as in Theorem \2. 11 then 



E° P o ;M (t)= Yl K0,u)E P o nu (t), (2.1) 

and if "K is transverse to P, 

E P ^(t)= \K0,u)\E Pnu (t). (2.2) 

ueL(P°,J{) 

Often the polytope is not full-dimensional. Suppose that s is any affine subspace. Its period p(s) 
is the smallest positive integer p for which p~ lr L d meets s. 

Corollary 2.3 ([1J Corollary 4.3]). Let P be a rational convex polytope and "K a hyperplane arrange- 
ment in s := aff P. Then Ep t j{(t) and E° pa ^(t) are quasipolynomials in t that satisfy the reciprocity 

lawE° po M (t) = (-l) dims E P ^(-t). Their period is a multiple ofp(s) and a divisor of the denomina- 
tor of(P,'K). Ift = mod p(s), the leading term of Ep^c(t) is (vol^-i^d p)t dims and its constant 
term is the number of regions of (P,!K); but ift ^ mod p(s), then Ep^(t) = E° po ^(t) = 0. ■ 



3. Magic squares, magic labellings of covering clutters, and equal line sums 

3.1. Sorts of magic. We pointed out in the introduction that the difference between weak and 
strong magic (or semimagic) squares lies in the fact that for the latter we require the entries to 
be distinct. We will therefore spend the beginning of this section studying the general setting of 
integer points in polytopes with distinct entries. 

We have a convex polytope P CM. d , spanning an affine subspace s. To ensure distinctness of the 
coordinates of a vector we avoid the hyperplanes of "K := J{[Kd] s , the complete-graph arrangement 
3-C[-?Q] intersected with s. Transversality of "K and P means that, first of all, s is not a subspace of 
any hyperplane Xj = x^ and, secondly, any flat of !H[-fQ] that meets P also meets P°. In many of 
the interesting special cases the latter condition is automatic. 

Suppose x is a point in M. d whose entries are all distinct. There is a unique permutation r of 
[d] such that x T \ < x T 2 < • ■ ■ < x TC i- We say x realizes r. We call a permutation a realizable in a 
subset A C s if there is a vector x £ A that realizes it. We are interested in realizability in P, but 
realizability in s is simpler. Fortunately, a permutation that is realizable in s is also realizable in 
P if P contains a positive multiple of 1 := (1, 1, . . . , 1), since every closed region of "K contains (1). 
This is the case when every form has equal positive weight. 

If x 6 R d satisfies x a \ < %ct2 < • • • < Xcrd, we say x and a are compatible. 

Theorem 3.1. Suppose P C M d is a closed, rational convex polytope transverse to 3~C[-?Q] and 
s := aff P. The number E° po j£^ d ys(t) of integer points in tP° with distinct entries is a quasipolyno- 

mial in positive integers t with leading term (vol P)t dims . Furthermore, (—l) dims E po ^^-j^ - 1) = 
Ep,w[K d ] s (t) := the number of pairs (x,a) consisting of an integer point x G tP and a compat- 
ible P -realizable permutation a of [d]. The constant term of E P ^ Kd ]s(t) equals the number of 
permutations of [d] that are realizable in P. 

Proof. The first statement is a direct consequence of Corollary 12 . 31 along with the observation that, 
by transversality, a region that intersects P must also intersect P°. For the second, the multiplicity 
m(x) of x equals the number of closed regions of iKf-fQ] 5 that contain x. The regions of JCf-fCd] 5 
correspond to certain regions of !K[K^], which correspond to permutations of [d\. Clearly, a closed 
region contains x if and only if its permutation is compatible with x. Thus, m{x) is the number 
of permutations that are both realizable in P and compatible with x. Now appeal to Corollary 
I2~31 ■ 
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Problem 3.2. The period and denominator present a puzzle. The denominator of the inside-out 
polytope (P,!K[Kd\) is obviously a multiple of the denominator of the standard polytope P. The 
first question is when the hyperplane arrangement "K [Kd] changes this latter denominator, and in 
what way. As for the period, if in particular P has integral vertices then Ep is a polynomial. What 
conditions on P ensure that Epw Kd ^ is also a polynomial? That it need not be is illustrated by 
the simple example of the line segment from (0, 1) to (1, 0) in M? and the hyperplane x\ = xi- 

Our desire to develop the ideas behind magic squares and graphs suggests two approaches to 
choosing P and s. The subspace, s, represents the existence of a magic sum. The polytope, P, 
represents the constraints on the entries in the value vector. The magic sum constraints may be 
pure equalities: 

(i) Set all linear forms equal to each other. (Homogeneous equations.) 
Or, they may be set all equal to a controlled constant: 

(ii) Set all linear forms equal to t. (Affine equations.) 

(Sometimes one wants additional equations; see the discussion of centrally symmetric squares in 
Examples 13.91 and 13. 201 ) Similarly, the constraints on the components of x may be two-sided 
bounds: 

(I) All variables Xi satisfy < Xi < t. (Cubical constraints.) 

Or, the constraints may be merely nonnegativity of the variables: 

(II) All variables xi > 0. (Nonnegativity.) 

We think the natural combinations are (i) with (I) and (ii) with (II) and that is how we develop 
the theory. 

3.2. Cubical magic. The cubical approach to magic squares counts them by the largest allowed 
value of the entry in a cell; if t is the parameter, the squares counted are those with entries 
< i/ij < t. Similarly, magic labellings of a bidirected graph are counted by the upper bound t — 1 
on the edge labels. 

In the general situation the magic subspace s is defined by homogeneous, rational linear equations 

/i(x) = / 2 (x) = ... = / m (x) = 0. (3.1) 

These are obviously equivalent to (i), if the forms in (|3.1[) are the differences of the forms of (i). 
The magic polytope is 

P := [0,l] d ns, 

whence the name "cubical". The hyperplane arrangement is !H := !K[ir^] ,s . In the example of magic 
squares, P is the set of all square matrices with real entries in [0, 1] such that all row, column, and 
diagonal sums are equal. For magic graphs P is the set of edge labellings by numbers in [0, 1] such 
that all node sums are equal. 

One has to be sure that P spans s. 

Lemma 3.3. // P is not contained within a coordinate hyperplane, it affinely spans the magic 
subspace s. 

Proof. By the hypothesis P contains points x 1 (not necessarily distinct) with < (x l )i < 1. The 
barycenter of these points lies in (0, l] d n s, hence in P. Therefore P° = s n (0, l] d , which clearly 
spans s. ■ 

If P is contained in a coordinate hyperplane, one should reformulate the problem with fewer 
variables. 

Now for the main theorem on cubical magic. The magic subspace s C R d is given by (|3.1|) : the 
polytope is P := [0, l] d (Is. For t = 1, 2, . . . let 
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B°(t) := the number of integer points x € s with distinct entries that satisfy < Xj < t, 
and let 

B(t) := the number of pairs (x,a) consisting of an integer point x G s that satisfies < X{ < t 
and a compatible s-realizable permutation a of [d]. 

Theorem 3.4 (Magic enumeration by bounds). Suppose P := [0, l] d fl s does not lie within a 
coordinate hyperplane. Then B° and B are quasipolynomials with leading term (vol P)t dims and with 
constant term B(0) equal to the number of permutations of [d] that are realizable in s. Furthermore, 

(_l)dim S £o(_ t ) _ B Q.y 

Proof. Theorem 13.11 shows that B = Ep-yi and B° = E° po ^ are reciprocal quasipolynomials. The 
constant term equals the number of permutations that are realizable in P. The same ones are 
realizable in s, because a permutation realized by x is also realized by ax for any positive real a. 
By choosing small enough a we can put ax into P or —P. If the latter, then ^1 — ax realizes the 
permutation and lies in P. ■ 

Most interesting is the case in which all the forms /j have weight zero. It is precisely then that "K 
and P are transverse, as one can see by comparing the subspace p| "K = (1) n s with the definition 
of transversality; moreover, then P spans s. 

Theorem 3.5. With s defined by forms of weight zero, and assuming a magic labelling exists, we 
have 

B°{t)= Yl A*(M%l)"n«(*) , 
«e£(P°,:K) 

B(t)= Y, IMM|£[o,i]<*n«(*) , 
«e£(P°,:K) 

where [i is the Mobius function of £ (P° , "K) . 

Proof. Apply Theorem 12.21 in s. A magic labelling exists if and only if s does not lie in any 
hyperplane Xj = Xk- ■ 

A flat u G £(!K) has the form v fl s where v is given by a series of equations of coordinates: 
Xjj = Xi 2 = ■ ■ ■ = Xi p , Xj x = Xj 2 = ■ ■ ■ = x j , etc.; that is, v corresponds to a partition ir of 
X. We can treat these equations as eliminating the variables Xj 2 , . . . , Xi , Xj 2 , . . . , Xj q , ... in favor 
of x^jXj!, .... With this substitution, v = R d ' for some d! < d and [0, l] d n v = [0, l] d ' . Then 
[0, l] d fl u is essentially [0, l] d fl s', where s' is s after identifying variables. Similar remarks apply 
to £(P°,IK), since its members are the nonvoid intersections with P° of flats of Ji. 

We describe some of the most interesting examples, concluding with the two best known. 

Example 3.6 (Lines of constant length). In a covering clutter (X, XL), suppose every line has the 
same number of points. Then the linear equations that express the existence of a magic sum take 
the form 

Xj = Xj for all Li, L2 € & ■ 

Theorems 13.41 and 13.51 both apply: B°(t) and B{t) are reciprocal quasipolynomials in the upper 
bound t — 1 or t, respectively, and so on. Examples include magic, semimagic, and pandiagonal 
magic squares, afline and projective planes, /c-nets, and magic hypercubes with or without diagonals 
of various kinds, as well as magic labelling of regular graphs. One has to ask about the existence of 
a magic labelling. There is no known general answer, but certainly there exist magic and semimagic 
squares of all orders n > 3 and pandiagonal magic squares of all orders n > 4, using the standard 
entries {1,2,..., n 2 } if n ^ 2 mod 4 (see % pp. 203-211]). 
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For a general covering clutter only Theorem 13,41 applies. If it has magic labellings at all, the 
problems of existence and characterization of realizable permutations come to the fore. An example 
of nonexistence is the Fano plane. Examples of existence are magic squares of size n > 3. Regarding 
characterization we propose a conjecture. We may assume X = [d]. Given a covering clutter, a 
magic permutation (with respect to the covering clutter) is a permutation of [d] that is realizable 
by a positive point x G P. Obviously, x can be chosen to be rational if it exists at all. In the 
cubical situation we are discussing here, all points in P° are positive, so magic permutations are 
identical to P-realizable and therefore to s-realizable permutations. A permutation a of [d] defines 
a reverse dominance order on the power set ^Qd]) by 

L =4 a L' if, when L and L' are written in decreasing order according to a, say L = 
{aji, ...,crji} where Ji > • • • > ji and L' = {aj[, aj' v } where j[ > ■ ■ ■ > j' v , then 
I < I' and ji <j[,...,ji <j[. 

This is a partial order on !P([(i]). 

Conjecture 3.7 (Magic permutations). A permutation a of [d] is realizable by a positive point in 
the magic subspace s of the covering clutter ([d], £) if and only if £ is an antichain in the reverse 
dominance order due to a. 

The archetype is magic squares. A magic permutation there is a permutation of the cells of the 
square that is obtained from some magic labelling by arranging the cells in increasing order. We 
have verified the conjecture for 3x3 magic and semimagic squares (Examples 13 . 1 1 1 and 13. 1 2]) . 

A permutation being a total order on [d], one could generalize to total preorders (in which the an- 
tisymmetric law is not required); we propose the corresponding conjecture, replacing "permutation" 
by "total preorder". 

The necessity for £ to be an antichain is obvious. It is the same with the extension to linear 
forms. We define the reverse dominance order due to a on the set of positive linear forms 

on R d by 

/ 4a f if, writing / = J2k=i a k x k and /' = Y?k=i a 'k x k, then 
d d 

^2a ak <^2 a 'ak for ever Y 3 = !> 2 > • • • > d. 

k=j k=j 

Conjecture 3.8 (Magic permutations for linear forms). Let /i, . . . , f m be positive forms on R d and 
s the subspace on which they are all zero. A permutation a of [d] is realizable by a positive point 
in s if and only if {f\, . . . , f m } is an antichain in the reverse dominance order of forms. 

Example 3.9 (Cubical symmetry). A cubically symmetric magic or semimagic square has the prop- 
erty that any two cells that lie opposite each other across the center have sum equal to t, the 
cell-value bound, and the center cell (if there is one) contains the value This definition is 
generalized from that of associated square in (symmetrical square in [2j), in which the entries 
are 1, 2, . . . , n 2 and a symmetrical pair sums to n 2 + 1. (Cubical symmetry contrasts with affine 
symmetry, which we shall treat shortly.) The novel feature is the additional linear restraints besides 
the magic sum conditions: these are j/y + y n+ \-i in+ \-j = t for all i, j G [m]. Translated into the 
language of ^-fractional vectors x = \y G \7L n , we require 

x ij ~\~ Xn+1— i,n+l— j = 1- (3-2) 

The effect is to reduce the magic subspace s to a smaller subspace s', but Theorems 13.11 and 13.51 
apply (with s' replacing s) if the line size is constant, due to the next lemma. In it we generalize 
symmetry to fairly arbitrary covering clutters on the point set [n] 2 . The definition is the same: we 
take the magic sum conditions and the symmetry equations ()3.2I) . The magic subspace s C R n is 
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defined by equality of line sums; the cubically symmetric magic subspace, s' , is the affine subspace 
of s in which (j32J) is valid; then P = [0, l]™ 2 n s' and "X = l K[K n 2] s> . 

Lemma 3.10. If f[n] 2 ,£) is a covering clutter of constant size such that s' is not contained in a 
hyperplane Xij = x/-i of 0i[K n a], and in particular if a cubically symmetric, strongly magic labelling 
exists, then P and are transverse. 

Proof. Since |l € PPif~]J{, every flat of "K intersects P°. ■ 

The hypothesis on s' is satisfied in the case of cubically symmetric magic squares of side n > 3 
because such squares are known to exist, using the values 1,2, ... ,n 2 if n ^ 2 mod 4; in fact, 
cubically symmetric pandiagonal squares exist if n > 4. See pp. 204-211]. 

Example 3.11 (Magic squares of order 3: cubical count). Define M°(t) for t = 1, 2, 3, . . . to be 
the number of 3 x 3 magic squares in which each cell value is less than t. From 0, the Ehrhart 
quasipolynomial is 

if t = 0,2,6,8 mod 12, 
if t = 1 mod 12, 
if* = 3, 11 mod 12, 
if * = 4, 10 mod 12, 
ift = 5,9 mod 12, 
if * = 7 mod 12. 



M°(*) 



t 3 


-16t 2 +76i- 


-96 _ 


(*- 


2)(t-6)(t-8) 




6 






6 


t 3 


-16t 2 +73t- 


-58 _ 


(t- 


l)(i 2 -15t+58) 




6 






6 


t 3 


-m 2 +73t- 


-102 _ 


(t 


-3)(t 2 -13i+34) 




6 






6 


t :i 


-m 2 +76t- 


112 _ 


(t 


-4)(t 2 -124+28) 




6 






6 


t 3 


-m 2 +73t- 


-90 _ 


(t- 


2)(t-5)(*-9) 




6 






6 


t :i 


-m 2 +73t- 


-70 _ 


(t- 


7)(t 2 -9i+10) 




6 






6 



The constant term of M c (t) = (— 1) 3 M°(— *) is the number of magic permutations, which is 16. 
These permutations are the rotations and reflections of the patterns 



(a) 



4 


9 


2 




3 


9 


2 


3 


5 


7 


(b) 


4 


5 


6 


8 


1 


6 




8 


1 


7 



In these diagrams the numbers are not cell values but rather permutation positions: the largest value 
is in the cell marked 9, the next largest in that marked 8, and so on. (To realize the permutations 
by magic squares, (a) can be left untouched but (b) needs numbers.) The general form of a magic 
square is, up to the eight symmetries and an additive constant on each value, 



-p 


a + P 


—a 


-{a -(3) 





a — (3 


a 


-(a + /3) 


P 



where a > [3 > and a ^ 2(3. If a > 2(3 we get the magic permutation (a); if a < 2(3 we get (b). 
This proves there are just 16 magic permutations. 

Example 3.12 (Semimagic squares of order 3: cubical count). Let S°(t), for * > 0, be the number 
of semimagic squares of order 3 in which every entry is less than *. [Sj has exact formulas. The 
constant term |Sc(0)| = 1296 = 6 4 equals the number of 3 x 3 semimagic permutations, that is, 
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magic permutations for semimagic squares of order 3, in agreement with Conjecture %TT\ We verified 
this by hand, finding all semimagic permutations of order 3, based on the fact that a normalized 
semimagic permutation is a linear extension of the partial ordering implied by the super normalized 
form of a semimagic square developed in [3]. 

3.3. Affine magic. The affine approach counts magic squares, and magic labellings in general, by 
the magic sum. In the general situation the magic subspace s is defined by a rational, nonhomoge- 
neous linear system 

/i(x) = / 2 (x) = --- = / m (x) = l (3.3) 
that we assume is consistent. The magic polytope P is the nonnegative part of s, that is, 

P := sHO, 

where O := K> , the nonnegative orthant, and the hyperplane arrangement is !K := ^[KdY ■ 

Affine magic is quite similar to cubical magic, but there is something new: one has to worry 
about boundedness of P. Obviously, P is bounded if the defining linear forms fi in (j3.3j) are positive 
and every variable appears in a form. If the latter fails, we are simply in the wrong dimension, so 
we make the overall assumption in this section that every variable appears in at least one form. (A 
covering clutter satisfies this automatically.) As in the cubical treatment, one must make sure that 
P affinely spans s (or else change s in the theorems to aff P) and one has to be concerned about 
transversality of P and !H. 

Lemma 3.13. If P is not contained within a coordinate hyperplane, it spans s. 

Proof. As with Lemma 13.31 ■ 

Let s C W d be the solution space of (|3.3j) . where the fi are rational linear forms, and let P := sHO. 
For t = 1,2,..., let 

A°(t) := the number of integer points x € tP with distinct positive entries, 
and let 

A(t) := the number of pairs (x,o~) consisting of a nonnegative integer point x £ tP and a 
compatible P°-realizable permutation a of [d\. 

When P is bounded and spans s, Theorem 13. II applies. 

Theorem 3.14 (Magic enumeration by line sums). Suppose that P is bounded and does not lie 
within a coordinate hyperplane. Then A° and A are quasipolynomials with leading term (vol P)t dims 
and with constant term A(0) equal to the number of permutations of [d] that are realizable in P° . 
Furthermore, (-l) dim s A° (-t) = A(t). 

Proof. A straightforward application of Theorem 13.11 ■ 

Problem 3.15. In Theorem 13.141 can realizability in P° be replaced by the weaker property of 
realizability in s? 

We want to know when P and "K are transverse so that the Mobius- function formulas of Theorem 
12.21 will apply, and also to allow a positive answer to the question of Problem 13.151 for which the 
combination of transversality and nonemptiness of P Pi (1) is necessary and sufficient. 

Lemma 3.16. P and "K are transverse if s is not contained in any hyperplane Xj = x^ and all 
forms fi have equal positive weight. 

Proof. We remarked that P is bounded because every coordinate appears in a form. We need to 
verify that for no u € is / P n u C dP. But P n u = P n v for some v € &(!K[K d ]), and 

dP is contained in the union of the coordinate hyperplanes since they determine the facets of P. 
If some v D P lies in a coordinate hyperplane, then the same is true of v = (1); but (1) fi P is in a 
coordinate hyperplane if and only if it is or {0}. The latter is impossible with forms as in ()3,3|) . 
If the forms have equal weig ht c > 0, then (1) n P = {c^l}. The converse is obvious. ■ 
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Theorem 3.17. With s defined by forms of constant positive weight, and assuming a magic labelling 
exists, we have 

A°(t)= M(M-EunO«(t) , 

M*)= \Ku,u)\E uno (t) , 

where fi is the Mobius function of £ (P° , J{) . 

Proof. Transversality holds by Lemma 13.161 since a magic labelling exists if and only if s does not 
lie in any hyperplane Xj = Xk- Apply Theorem 12.21 in s. ■ 

Some interesting examples are the affine versions of the cubical examples we already mentioned. 
In many of them the affine intersection poset is combinatorially equivalent to the cubical intersection 
poset. Let P c be the cubical polytope s n (0, l) d ; P is still the affine polytope sf)0. 

Lemma 3.18. Suppose that s is defined by positive forms such that each variable xa has some 
form in which its coefficient is at least 1. Assume also that a magic labelling exists. Then the 
affine and cubical intersection posets are naturally isomorphic: L{P° , ( K\K ( j\) = L(P° ,"H{Kd\) by 
uC\P° HiinP° for u e L(J{[K d }) such that unP c ^ 0. 

Proof. The assumption on the coefficients implies that P = P c n s. M 

Example 3.19 (Lines of constant length, cf. Example 13. 6|) . The linear equations that express the 
existence of a magic sum t take the form 

Xj = t for all L € £ . 

Theorems 13 . 1 41 and 13 . 1 7l and Lemma 13.181 all apply as long as a magic labelling exists; thus, A°(t) 
and A(t) are reciprocal quasipolynomials in the magic sum t, etc. 

For general covering clutters Theorem 13.141 applies, showing that A(t) and A°(t) are reciprocal 
quasipolynomials in t — provided that magic labellings exist at all. With that assumption, Lemma 
13.181 also applies. The affine magic subspace, call it s±, defined by X^ei xi = 1 for all L E £, is 
an affine subspace of the homogeneous magic subspace, call it so, and so is the linear subspace 
generated by s\. This means that the permutations realizable in sq and s± are the same, so all 
our comments on magic permutations with respect to a covering clutter, in the context of cubical 
counting, apply as well to affine enumeration, except that we do not in general know that all 
s-realizable permutations are P-realizable. 

Example 3.20 (Affine symmetry; cf. Example 13. 9|) . An affinely symmetric magic or semimagic 
square has the property that the average value of any two cells that lie opposite each other across 
the center equals the average cell value, t/n, and the center cell (if there is one) contains the 
value t/n. (Of course, one cannot expect such squares to exist unless t = (mod n).) This 
definition is another generalization of that of associated square. The additional linear restraints 
are yij + y n +i~i,n+i-j = 2t/n for all i,j£ [m]. Translated into the language of ^--fractional vectors 
x = \y € jZ n2 , we require 

2 

n 

The effect is to reduce the magic subspace s to a smaller subspace s'. With the extra hypothesis that 
all lines have n points, Theorems 13.11 and 13.171 apply, with s' replacing s, by the following lemma, 
in which we generalize symmetry to covering clutters on X = [n] 2 , with magic sum conditions and 

the symmetry equations (]3.4[) . The magic subspace s C M n is defined by equality of line sums; 

li 



the affinely symmetric magic subspace, s', is the affine subspace of s in which (|3,2j) is valid; then 
P = R| a n s' and 5f = l K[K n 2] s ' . 

Lemma 3.21. If {[n] 2 , £) zs a covering clutter of size n such that s' is not contained in a hyperplane 
Xij = Xki of Ji[K n 2], and in particular if an affinely symmetric magic labelling exists, then P and 
"K are transverse. 

Proof. Since ~1 G P n f| "K, every fiat of "K intersects P°. ■ 
Example 3.22 (Magic squares of order 3: affine count). Let M°(t), for t > 0, be the number of 



Kit) 



2t 2 


-324+144 




9 


2t 2 


-324+78 




9 


2t 2 


-324+120 




9 


2t 2 


-324+126 




9 


24 2 


-324+96 




9 


2t 2 


-324+102 



j5J we find 


that 








-lit 2 


- m 


+ 72) 


if t = 





mod 18, 


lit- 


3)(t- 


-13) 


if t = 


3 


mod 18, 


-lit- 


-6)(i 


-10) 


if t = 


6 


mod 18, 


-w- 


-m 


-9) 


if t = 


9 


mod 18, 


lit- 


m- 


-12) 


if t = 


12 


mod 18, 


-lit 2 


- m 


+ 51) 


if t = 


15 


mod 18, 



if t ^ mod 3. 

Xin ^H] has another way to find the generating function of M°(t) (with the minor difference 
that he allows zero entries), by applying MacMahon's partition calculus. 

The constant term M a (0) = 16 is the number of magic permutations; this is the same number 
as with cubically counted magic squares, M c (0) in Example 13.111 

Example 3.23 (Semimagic squares of order 3: affine count). Complete results for this example are 



m 



0. 



4. Latin squares join in magically 

The general picture that encompasses latin squares is that of a covering clutter (X, £), as in 
Section |31 with an integer labelling x : X —> Z subject to the requirement that 

x{e) ^ x(f) if e and / lie in a line. 

This is a latin labelling of (X, £). The graph of forbidden equalities is therefore 

IY := (J K L , 

Kl being the complete graph with L as node set. Every graph is equal to T^, for some choice of 
covering clutter. An orientation of T^ is acyclic if it has no directed cycles. An orientation and a 
node c-coloring x : V — > [c] are compatible if Xj > Xi whenever there is a T^-edge oriented from i 

to j mi. 

A crucial decision is how to restrict the symbols of the latin labelling. One may simply specify 
the number of symbols allowed, say k, and an arbitrary symbol set, let us say [k]. Then the number 
of latin labellings equals the chromatic polynomial xr^ {k), so we would have merely an application 
of graph coloring. In this article, taking a leaf from the guidebook of magic squares, we focus on 
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approaches which give rise to magilatin squares: we impose a summation condition on the lines, 
which may be either homogeneous: 

(i) Set all line sums equal to each other, 
or affine: 

(ii) Set all line sums equal to t. 

This tactic is not so strange as it may appear. A partial magilatin square with homogeneous line- 
sum requirements and symbols restricted to the interval [1, n] is just a latin square with symbol set 
[n]. So is a partial latin square with affine line-sum restrictions, line sum t = C^ 1 ), and positive 
symbols. We are assuming that the affine constraint (ii) is supplemented by a positivity assumption 
and that the homogeneous constraint (i) is supplemented by the requirement that the symbols be 
drawn from the set [k] for some k. To handle partial latin rectangles requires a generalization 
of (i) to multiple covering clutters; see Section 14.11 In every case the hyperplane arrangement is 
!K := !K[r^] s , where s is the subspace of R determined by the appropriate line sum conditions. 

4.1. Cubical latinity. The cubical approach to latin labellings counts them by the largest allowed 
value; if t is the parameter, the labellings counted are those with < x\ < t. We assume a multiple 
covering clutter, (X; &i, . . . , that is, there are k > 1 covering clutters in our 

examples k will be 1 for squares and 2 for rectangles. The labelling must be such that the lines in 
each class have equal sums; but the sums in different classes are independent of each other. For a 
magilatin labelling the numbers in each line must be distinct. To count the labellings we take the 
subspace s of M. d in which all line sums within each covering clutter are equal. The polytope 
is P := [0, l] d n s. This is as in Section 13.21 so Lemma 13.31 applies to assure that P spans s. For 
t = 1,2,..., let 

L°(t) := the number of latin labellings x with equal line sums within each covering clutter and 
with entries that satisfy < X{ < t, 

and let 

L c (t) := the number of pairs consisting of a latin labelling x, with equal line sums within each 
covering clutter and with < x% < t, and a compatible acyclic orientation of T^,. 

Theorem 4.1 (Magilatin enumeration by upper bound). Suppose P := [0, 1] Ds does not lie within 
a coordinate hyperplane. Then L° and L c are quasipolynomials with leading term (volP)t dims and 
with constant term L c (0) equal to the number of acyclic orientations ofT& that are realizable in s. 
Furthermore, {-l) dims L°{-t) = L c (t). 

Proof. Theorem 12.11 shows that L c = Ep^ and L° = E po ^ are reciprocal quasipolynomials. The 
remainder of the proof is as in Theorem 13.41 ■ 

Theorem 4.2. If all lines have the same size and a latin labelling exists, then we have all the 
formulas of Theorem \3.5\ with L c in place of B . 

Proof. Apply Theorem 12.21 ■ 

Example 4.3 (Small latin shapes, counted cubically). We calculated two very small magilatin ex- 
amples: a square and a rectangle. 

In the 2x2 magilatin square it must be that x\\ = X22 7^ £12 = ^21- With cubical constraints, 
then, 

L c °(t) = (t-l)(t-2). 

For comparison, the number of squares without the magilatin distinctness requirement is t 2 . 

The 2x3 latin rectangle is much more complicated. Calculation with Maple shows that the 
denominator of P is 6 and that of (P, "K) is 12. This permits us to calculate the quasipolynomial 
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(by actual count of rectangles up to t = 60). Its period turns out to be 4. 

> t3-12^+41 f -30 = (t-l)ft-5)( t -6) if t = ! mod4) 
L° c (t) = \ * 3 -12f 2 +41f-42 = (t-2)(*-3)(*-7) if f _ 3 mod ^ 
^-12f+44t-48 = (t-2)( t -4)(«-6) if t ig eyen . 

For comparison, the number of rectangles without the distinctness requirement is 

f f 3_ 3f 2 +3f _ 1 _ (£-1)3 . f 4 ig odd ^ 



4 ~~ 4 

if t is even. 



f 3 -3f 2 +6f-4 _ (f-l)(t 2 -2t+4) 



4 — 4 

Problem 4.4. We do not know why the coefficients alternate in sign, what causes the smallness of 
the differences among the constituents, nor what makes the even constituents have a smaller period 
than the odd constituents. 

Example 4.5 (Magilatin 3x3 squares, counted cubically). See jS| for the complete solution. 

4.2. Affine latinity. For affine counting of latin labellings we take a covering clutter (X,Ju). s is 
the subspace in which all line sums equal 1, and P = s n O, as in affine magic (Section I3.3[) . As 
there, we can apply Lemma 13.131 to conclude that P spans s in all interesting cases. One advantage 
over magic is that, because we consider only line sums and not general linear forms, P is certain 
to be bounded. For t = 1, 2, . . ., let 

L°(t) := the number of latin labellings of (A, £) with positive entries and all line sums equal to 
t, 

and let 

La,(t) := the number of pairs consisting of a nonnegative latin labelling with all line sums equal 
to t and a compatible acyclic orientation of V n that are realizable in P°. 

Lemma 4.6. P and J£ are transverse if a positive latin labelling exists. 

Proof. As with Lemma f3.16l The existence of a positive latin labelling implies that s is not contained 
in any hyperplane Xj = Xk for which a and ej are collinear and, for every v 6 £(Jf[r^]), v fl P is 
not contained in a coordinate hyperplane. ■ 

Theorem 4.7 (Magilatin enumeration by magic line sum). Suppose that a positive latin labelling 
exists for some t. Then L° and L a are quasipolynomials with leading term (vol P)t dims and with 
constant term L a (0) equal to the number of acyclic orientations of that are realizable in P. 
Furthermore, (-l) dims L°(-i) = L a (t). 

Proof. By Theorem 13.11 adapted to the latin nonequalities, together with Lemma 14.61 to ensure 
transversality so that realizability in P is equivalent to realizability in P°. ■ 

If every line has the same size, then the acyclic orientations that are realizable in P° are the 
same as those realizable in s, because then (1) intersects P°; see the discussion following Theorem 

Theorem 4.8. Assuming a positive latin labelling exists and all lines have the same size, we have 
the formulas of Theorem \3. 1 7| with L a in place of A. 

Proof. Transversality holds by Lemma 14.61 Apply Theorem 12. 21 in s. ■ 

Example 4.9 (Small magilatin squares, counted affinely). It is easy to see that in a 2 x 2 magilatin 
square with affine constraints, 

t — 1 if t is odd, 
t — 2 if t is even. 



Kit) 



The period of 2 equals the denominator of (P, "K). The number of positive squares without the 
magilatin distinctness requirement is t — 1. 
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Example 4.10 (Magilatin 3x3 squares, counted amnely). For the number of these squares see 

5. Generalized exclusions 

We concentrated our treatment on exclusions of magic and latin type: that is, where all values, 
or all collinear values, are unequal. Other kinds of exclusion are possible. We want to mention the 
very natural complementation restrictions. In cubical enumeration we call Xi and Xj complementary 
if Xi + xj = t. If we forbid certain pairs of values to be complementary, we pass from graphs to 
signed graphs, since the rule Xi + Xj ^ t corresponds to a negative edge —ij. An inequality Xi ^ Xj 
corresponds to a positive edge, thus ordinary edges are positive. The exact application of signed 
graphs involves translation and halving of the centrally symmetric polytope [—1, l] d to [0, l] d , along 
with corresponding translation of the signed-graphic hyperplane arrangement, as explained in [IJ 
Section 5]. Because signed-graphic hyperplanes, as translated, give half-integral vertices, we expect 
a counting quasipolynomial with nonequalities and noncomplementarities given by a signed graph 
to have twice the period of a counting quasipolynomial pertaining to a similar unsigned graph 
of inequalities. This is necessarily vague; we intend only a suggestion for research that we invite 
readers to explore. 
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